Casimir interaction between a dielectric nanosphere and a metallic plane 
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We study the Casimir interaction between a dielectric nanosphere and a metallic plane, using the 
multiple scattering theory. Exact results are obtained with the dielectric described by a Sellmeier 
model and the metal by a Drude model. Asymptotic forms are discussed for small spheres, large 
or small distances. The well-known Casimir-Polder formula is recovered at the limit of vanishingly 
small spheres, while an expression better behaved at small distances is found for any finite value of 
the radius. The exact results are of particular interest for the study of quantum states of nanospheres 
in the vicinity of surfaces. 



PACS numbers: 

I. INTRODUCTION 

The Casimir effect, due to the scattering of quantum 
fluctuations of the electromagnetic vacuum is the 
dominant interaction between neutral bodies at distances 
large compared to atomic scales 0, Q (more references 
in [j]). For this reason, it has a strong impact in var- 
ious important domains, such as atomic and molecular 
physics, condensed matter and surface physics, chemical 
and biological physics, micro- and nano-technology 

In this paper, we will consider the case of dielectric 
nanospheres in the vicinity of a metallic surface. This 
study is part of a discussion of the intriguing phenomenon 
of small heating of ultra-cold neutrons (UCN) in traps @- 

which could be explained by the interaction between 
UCN and nanospheres levitated in the quantum states 
created by the interaction of nanospheres with surfaces 
0. In order to characterize this phenomenon and com- 
pute the properties of the quantum states, one needs to 
have a detailed and careful treatment of the interaction 
potential. In particular, as shown below, the commonly 
used Casimir-Polder formula is not sufiicient to this 
purpose. 

Below we first recall how the scattering formula [Tl| 
can be applied to the study of a sphere of radius R at 
a distance L of closest approach to the plane. We then 
give numerical evaluations and graphical plots of the in- 
teraction energy. The Casimir-Polder formula is recov- 
ered when the radius of the nanosphere is smaller than 
all other length scales. The short and long distance limits 
are then found to show subtle interplays with the limit 
of small radius. 

Important consequences of these results are obtained 
for the behavior of the interaction at small distances. 
While the Casimir-Polder energy scales as L^^ at small 
values of L, the full expression is found to be better be- 
haved for any finite value of R, which leads to a regular 
solution for the quantum states. 



II. SCATTERING FORMALISM 

In this paper, we will not consider the effect of thermal 
fluctuations on the interaction (they are expected to be 
small at the not so large separations considered here). 
We thus start from the scattering formula for the Casimir 
energy at zero temperature [ill ] 

E^hj ^IndetP, V=[I-M) (1) 

M=nse-^^npe-^^ , C = L + R 

The Casimir energy is written in terms of reflection op- 
erators TZs and TZp which describe the diffraction by the 
sphere and the plate. These operators are evaluated with 
reference points at the sphere center and at its projection 
on the plane, respectively. The operator e~'^^ accounts 
for one-way propagation along the distance C = L + R 
separating these two points. The operator M thus repre- 
sents one round-trip propagation inside the cavity formed 
by the two surfaces. All quantities are written at imagi- 
nary frequencies lj = after a Wick rotation. 

The scattering formula ^ provides a compact way of 
taking the multiple scatterings between the interacting 
bodies into account. It can be considered as generalizing 
the Dzyaloshinskii-Lifshitz-Pitaevskii formula [l2] to ar- 
bitrary scattering properties of the two bodies. It can be 
applied in various geometries and has in particular been 
recently used for calculating the Casimir interaction be- 
tween a metallic sphere and a metallic plane [l3l - [T6l | (see 
also [13, [13 )• III the following, we use the same tech- 
niques and notations as in [l5|, [lQ| and apply them to the 
case of a dielectric nanosphere above a metallic plane. 

The reflection on the plane is conveniently written by 
using a plane- wave basis |k, ±,p) where k is the wavevec- 
tor component parallel to the plane xy of the metal- 
lic surface, p = TE,TM the polarization, +/— the up- 
wards/downwards propagation direction. This basis is 
well adapted to the description of the propagation op- 



2 



erator e" 



-KC 



since /C is thus diagonal with elements 



AC = l<? + k"^ representing the (Wick rotated) wave- 
vector z-componcnt for the imaginary frequency ^. The 
reflection operator p preserves all plane wave quantum 
numbers but the propagation direction, and its elements 
are given by the standard Fresnel specular reflection am- 
plitudes for an homogenous medium. 

Then the reflection on the sphere is more easily writ- 
ten by using the multipole basis |£mP), with -f- 1) 
and m denoting the angular momentum eigenvalues (with 
I = 1,2,..., m = ...,£) and P ~ E,M representing 
electric and magnetic multipoles. The reflection opera- 
tor TZs has its elements given by the standard Mie scat- 
tering amplitudes. Thanks to rotational symmetry, the 
operator Ai commutes with the angular momentum op- 
erator Jz. Hence Ai is block diagonal, and each block 
^("i) (corresponding to a given subspace m) yields an 
independent contribution to the Casimir energy 



fc poo ' 
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The primed sum is a sum over positive integers with the 
term m = counted for its half. We have organized 
Ai in terms of block-matrices built up on electric and 
magnetic contributions. The corresponding matrix ele- 
ments are given and discussed in |16| . 

We now apply these formulas to the case of a dielec- 
tric nanosphcrc above a metallic plane. In particular, the 
plots shown below will be calculated for the case of inter- 
est for UCN studies @, namely a diamond nanosphere 
above a copper plane. We model copper dielectric re- 
sponse with a Drude model 




FIG. 1: Relative permittivity e(i^) for copper (Drude model; 
red line) and diamond (Sellmeier model; black dashed). 



damping is disregarded here since it do not play any sig- 
nificant role. 

The two permittivities are shown on FiglTJ where the 
black dashed curve is for diamond, and the plain red 
curve for copper. For copper, s{i^) is very large for 
^ much smaller than wp, which means that the metal 
tends to become a very good reflector. For diamond, 
e(i^) tends to its static value e(0) = 1 -I- i?i for <^ much 
smaller than wi. For larger frequencies in contrast, the 
dielectric properties become poorer for diamond as well 
as copper. As and wi have similar values, we thus ex- 
pect a transition to take plac e between van der Waals and 
Casimir-Polder regimes when the distance between 
two objects is of the order of Ap or Ai. The simple models 
(|3]|4|) are sufficient for the purpose of the present work. 
They could easily be improved to take into account in- 
terband transitions for copper, multiple components and 
damping in the Sellmeier model for diamond. 



eitO = 1 + 



(3) 



This is written at imaginary frequencies ut = with 
the squared plasma frequency proportional to the density 
of conduction electrons in the metal, and 7 the damping 
rate which measures the relaxation of these electrons. For 
explicit calculations and plots, we will use the relations 
Ljp = l-KcjXp with the plasma wavelength Ap = 136nm 
and 7 = O.OOSSwp. As 7 is small when compared to ujp 
for a good metal such as copper, its influence is small 
at the not too large distances considered in the present 
study (see [l^ for more details and references). The 
diamond dielectric response is described by a Sellmeier 
model 



£(*e) = i + E 



one 



(4) 



For diamond, a good enough description will be obtained 
with a single component in this formula with B\ = 4.91 
and oJi = l-KcjXx with the wavelength Ai — 106nm. The 



III. NUMERICAL EVALUATIONS 

For evaluating the determinant in (IT][5]), one needs to 
truncate the vector space at some maximum value ^max of 
angular momentum |13l - ll6l | . A qualitative understanding 
of the associated effects may be obtained from the local- 
ization principle [20j : the value of £niax required for a 
given accuracy level is expected to scale with the size 
parameter ^i? (where ^ = ^/c) which captures the de- 
pendence of scattering amplitudes on the sphere radius. 
Meanwhile, the frequencies giving the main contribution 
to the Casimir energy scale as ^ '--^ 1/L. As a conse- 
quence, the required ^max scales as RjL for intermediate 
and short separation distances. In the present paper, we 
will be interested in nanospheres i? < 20 nm so that cal- 
culating with £niax = 100 will be sufficient for a good 
accuracy for the results discussed below. 

The numerical results are shown in Figl2] for differ- 
ent values of R (2, 5, 10, and 20 nm). We plot the 
absolute value \E\ of the Casimir energy (E < 0) with 
respect to the distance L, from 1 to 500 nm. We see 
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FIG. 2: Absolute value of the Casimir energy E (measured 
in eV) with respect to the distance L for nanospheres of radii 
R = 2,5, 10, 20 nm. Logarithmic scales are used on both axis. 



that two regimes appear, which are reminiscent of non 
retarded van dcr Waals regime at short distances [2l| . 
and retarded Casimir-Polder regime at large distances 
, with different dependences upon the distance L (see 
the slopes of the curves) and the radius R (see the ver- 
tical spacing between the curves). Our results do not 
coincide with these well known limits because they take 
into account higher order multipole contributions to the 
scattering upon the sphere [l^. Of course, the Casimir- 
Polder expression will be recovered for small values of the 
radius, as shown in the next section. 



IV. LIMIT OF SMALL NANOSPHERES 



In this section, we consider the limit of a punctual 
sphere, when the radius R is smaller than any other 
length scale (in loose terms, we take the limit i? — 0). 
We show that the Casimir-Polder expressions are recov- 
ered, as expected. As the derivations presented here are 
applied to dielectric nanospheres, we will use the fact 
that the permittivity for the nanosphere remains finite 
at all frequencies. We also keep the same description of 
the dielectric response for the small values of the radius 
considered here. Note that for metallic nanospheres in 
contrast, it would be necessary to take into account the 
confinement effect for conduction electrons [l^, [13] . 

We start from the formulas ((TJ^l) giving the Casimir 
energy for a sphere of radius i? at a distance L of closest 
approach to a plane. Reflection on the plane is described 
by Fresnel amplitudes tte and ttm while scattering on 
the sphere is described by Mie amplitudes and (de- 
fined as in [1^). The following expressions, valid at low 
values of the parameter i?^, are sufficient for calculating 



the energy at the limit R ^ 
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ie + e + 1 {2e + iy.\{2e - ly.i 

(e- 1) [Rn 

^ ' {2e + 3y.\{2£ + iy.\ 



(5) 



As the dimensionless number R^ is much smaller than 
unity (for ~ 1; see the discussion in the previous sec- 
tion) and e remains finite at all frequencies, it follows that 
all these amplitudes are small and, simultaneously that 
the amplitude ai dominates all other Mie amplitudes. 

The calculation of the energy is therefore much simpler 
than in the general case. As a first simplification, one 
may indeed replace the non linear expression (0) by a 
linearized one (perturbative approximation) 
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Then, one may keep only the contributions to this sum 
which are proportional to the amplitude ai (electric dipo- 
lar approximation) 
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We also rewrite ai in terms of a dynamical electric polar- 
izability a(^) defined for the small nanosphere {aR^ is a 
reduced polarizability having the dimension of a volume; 
the SI polarizability is e^aR^) 



ai 
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(8) 



Collecting the results, we finally recover the full Casimir- 
Polder formula [l^l as written in [2^, \2l\ 
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We have used the fact that tte (iC) < and ttm (iC) > 0. 

We repeat at this point that the formula (jH]) has 
been obtained after two simplifications corresponding to 
the perturbative approximation and electric dipolar ap- 
proximation. As the Casimir-Polder interaction between 
atoms (loj . it contains in particular the limits of non 
retarded van der Waals and retarded Casimir-Polder ex- 
pressions [2^, 113] ■ The energy scales in both cases as R^ 
that is also the volume of the sphere. This result means 
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that the nanosphere behaves at the hmit i? as a big 
atom with an electric polarizabihty aR^. As we will see 
in the next section, this simple behavior does not remain 
true for arbitrary values of the radius R. 



V. ASYMPTOTIC BEHAVIOR AT SHORT AND 
LONG DISTANCES 

In this section, we discuss the asymptotic behaviors of 
the Casimir energy E at short and long distances. 

As a first step to this aim, we write a Casimir-Polder 
formula Eqp deduced from Ei in the limit L ^ Ap,Ai 
where copper may be considered as perfectly reflecting 
and diamond as having a constant electric polarizabihty 
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We proceed similarly with the van der Waals prediction 
i?vdw deduced from Ei for L <C Ap, Ai 



vdW 



C3 



3L3 



(11) 



16(V^Ap + VI -aoAi) 
The values i?cp ^^nd E^dw a-re equal at a crossing length 
3 {V2\p + VI -aoAi) 
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(12) 



that is approximately 39nm with the values correspond- 
ing to copper and diamond. 

As expected from already presented qualitative argu- 
ments, the exact Casimir expression E (given by eq[T]) is 
well approximated by Ecp when R is the smallest length 
scale and L the largest one 



E ~ Ecp , R<^L^<^L 



(13) 



Meanwhile, the exact Casimir expression E is well ap- 
proximated by Evdw when R is the smallest length scale 
and L smaller than L» 
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R-^L<^L^ 



(14) 



Of course, there exist a variety of behaviors when the two 
latter conditions are not met. 

In order to explore this variety, we plot on Fig|3] the 
logarithmic slope (log- log derivative) of the energy \E\ 
versus distance L 
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The parameter v would be a constant if the energy \E\ 
obeyed a power law dependence l/i" (for example j/ = 4 




10 
L [nm] 



1000 



FIG. 3: Logarithmic slope v as a, function of the distance 
L. Plain curves represent the nanosphere case with different 
values for the radius R, dashed curve represents the atomic 
limit, that is a nanosphere with _R — )■ 0. 



for f^cp or = 3 for i?vdw)- In the general case, v 
depends on L and can thus be understood as describing 
a 'local' power law in the vicinity of L. 

We see on Fig. [3]that v tends to the expected value 4 at 
large distances. At small distances in contrast, the van 
der Waals value = 3 is never a good approximation, 
which can be understood by inspecting the conditions 
for to be true. For any finite value of the radius 
R, we have indeed to cross the conditions L ^ R when 
the distance is decreased and (fT4)) can no longer be valid 
after this crossing. 

Another facet of the same problem becomes apparent 
when looking at the dependence of E versus R. We plot 
on FigH] a logarithmic slope calculated as in (fT3)) 
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The parameter fi would be a constant if obeyed a 
simple power law dependence R^. In particular the vol- 
umetric value /i = 3 is obtained for both Eqp and i?vdw- 
We see on Fig. 2] that fj, approaches this value at the limit 
of large distances, but departs from it everywhere else, 
anew indicating that E^dw is not a good approximation. 

We now give an improved version of the van der Waals 
formula ([TT|) which explains some of the features of the 
exact energy E. As i?vdw has been demonstrated above 
in the limit of a punctual sphere i? — > 0, wc may improve 
it for a finite size of the sphere through a pairwisc sum- 
mation over the volume. We thus obtain the Hamaker 
expression [28| 
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FIG. 4: Logarithmic slope ^ for variation with the radius R, 
as a function of the distance L. Plain curves represent the 
nanosphere case with different values for the radius R, the 
atom case being in comparison a constant = ?>. 

For completeness, we did proceed similarly with the 
Casimir-Polder formula -Bcp- 

These results allow one to understand the behaviors 
apparent on Figs[21lll Let us again consider that wc start 
from small nanospheres i? ^ at large distances L ^ 
L*. Using the expression i?cPj one obtains v = A and 
/X = 3. When the distance L is decreased, we cross two 
transitions L L^: and L ^ R and end up with the 
formula E^dw for which we get 1^=1 and fj. = 1. This 
line of reasoning reproduces the global variations seen on 

Figsum 
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FIG. 5: Ratios of the exact energy E to the expressions Ecp 
(plain blue curve) and -Evdw (dashed green curve). The two 
curves are calculated for R =10nm. 

In order to assess the quality of the estimations ([TT]) , 
we now plot on Figl5]the ratios E/Ecp and E/E^^w- As 
expected, we find that Eqp tends to reproduce the result 



E of the full numerical computation at large distances. 
We also see that i?vdw obeys the same power law than 
E at small distances (ratio tending to a constant value) , 
but fails to predict the correct magnitude (the limit of the 
ratio is not 1). This feature can be understood through 
a close inspection of the case where L is the smallest of 
all length scales. 

When L <^ R, wc can use the proximity force approx- 
imation and express the plane-sphere result in terms of 
the plane-plane one [2^ . We are thus left with the evalu- 
ation of the Casimir effect between copper and diamond 
plane plates. As L ^ L*, it is then possible to use the 
method designed in [s^l to find 

R 

EpFA^-TTc'sj, L^R,L, (18) 

This expression shows the same functional dependen- 
cies than the Hamaker expression E^dw ^ —nc^R/L 
while giving a different proportionality constant. This 
difference is due to the fact that P7)) has been ob- 
tained through a pairwise integration of van der Waals 
forces whereas has been calculated by taking into 
account the multiple interferences occuring in the Fabry- 
Perot cavity • The numerical values thus obtained 
(cg = 0.84 C3 with the numbers corresponding to dia- 
mond and copper) explain the behavior seen on Fig|SJ 

VI. CONCLUSIONS 

In this paper, we have presented an exact calculation of 
the Casimir interaction between a dielectric nanosphere 
and a metallic plane, using the multiple scattering for- 
malism as developed recently for the plane-sphere geom- 
etry [l3l - fl6j . In order to discuss qualitatively the results 
obtained in this manner, we have also investigated the 
limits of a punctual sphere as well as the asymptotic be- 
haviors at short and long distances. 

This study has important applications for discussing 
the intriguing phenomenon of heating of ultra-cold neu- 
trons in traps p-Q- This heating could be explained by 
the interaction between UCN and nanospheres levitated 
in the quantum states created by their interaction with 
surfaces In order to characterize quantitatively this 
phenomenon, a detailed knowledge of the interaction po- 
tential is required. In particular, the small-distance be- 
havior of the Casimir energy plays a crucial role in the 
determination of the quantum states obtained by solv- 
ing the Schrodinger equation for the wavefunction of the 
nanosphere in this potential. 

The commonly used Casimir-Polder formula, which 
also corresponds to the limit of our calculations for a 
vanishingly small radius R, leads to significant difficul- 
ties since it predicts a power law \E\ cx R^/L^ in the 
vicinity of the surface and thus leads to a ill-behaved 
Schrodinger problem. The exact solution presented in 
the present paper for a finite value of the radius R pre- 
dicts a smoother power law \E\ oc R/ L in the vicinity of 
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the surface and thus leads to a regular solution for the 
Schrodinger equation [91] ■ 

The plate roughness has been disregarded here and it 
is treated in a phenomenological manner in Q . It would 
be interesting to analyze its effect by using techniques 
already developed for treating the scattering on rough 
surfaces |3ll - [33| . 

It would also be worth investigating the same problem 
for the interaction between an atom and a plane. In 
analogy with the discussion of the present paper, taking 
into account the higher order multipoles and multiple 
interferences could lead to an expression of the energy 
more regular than with the commonly used electric dipole 



approximation. 
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